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ABSTRACT 

Recent measurements of Si IV/C IV ratios in the high-redshift Lya forest (Songaila 
& Cowie, AJ, 112, 335 [1996a]; Savaglio et ah, A&A, in press [1997]) have opened 
a new window on chemical enrichment and the first generations of stars. However, 
the derivation of accurate Si/C abundances requires reliable ionization corrections, 
which are strongly dependent on the spectral shape of the metagalactic ionizing 
background and on the “local effects” of hot stars in nearby galaxies. Recent models 
have assumed power-law quasar ionizing backgrounds plus a decrement at 4 Ryd to 
account for He H attenuation in intervening clouds. However, we show that realistic 
ionizing backgrounds based on cosmological radiative transfer models produce more 
complex ionizing spectra between 1-5 Ryd that are critical to interpreting ions of Si 
and C. We also make a preliminary investigation of the effects of He H ionization front 
non-overlap. Because the attenuation and re-emission by intervening clouds enhance 
Si IV relative to C IV, the observed high Si IV/C IV ratios do not require an unrealistic 
Si overproduction [Si/C > 3 (Si/C)©]. If the ionizing spectrum is dominated by “local 
effects” from massive stars, even larger Si IV/C IV ratios are possible. However, 
unless stellar radiation dominates quasars by more than a factor of 10, we confirm the 
evidence for some Si overproduction by massive stars; values Si/C « 2 (Si/C)© fit the 
measurements better than solar abundances. Ultimately, an adequate interpretation 
of the ratios of C IV, Si IV, and C H may require hot, collisionally ionized gas in a 
multiphase medium. 

Subject headings: galaxies: quasars: absorption lines- galaxies: intergalactic 
medium-galaxies: evolution 


1. Introduction 

The detection of metals associated with the “high” column density [N(H I) > 10^^'^ cm“^] 
Lya forest of absorbers in quasar spectra (Cowie et al. 1995; Tytler et al. 1995) provides new 
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challenges to the understanding of the nature and origin of these absorbers. The contamination 
by metals must now be accounted for in models for cloud evolution. It is also possible to perform 
a detailed analysis of the thermodynamic and ionization state of these absorbers, which have 
previously been probed only by their hydrogen Lya absorption and the integrated effect of their 
He II Lya absorption (Giroux, Fardal, & Shull 1995). If the clouds are photoionized, the detection 
of C II, C IV and Si IV may constrain the shape of the metagalactic radiation spectrum and 
indicate the epoch of dramatic changes in the ionizing radiation field at z > 3 (Songaila & Cowie 
1996a; Savaglio et al. 1997). 

The multiple C IV lines associated with a single H I absorption complex (Songaila & Cowie 
1996a) belies a simple picture of a homogeneous cloud. These absorbers, like Lyman Limit 
systems, may better be interpreted with multiphase models (Giroux, Sutherland, & Shull 1994; 
Petitjean, Riediger, & Rauch 1996; Haehnelt, Rauch, & Steinmetz 1997; Hellsten et al. 1997) that 
may include hot, collisionally ionized gas. We defer such interpretation to a later study, focussing 
here on the tradition of analyzing metal-line ratios with single-phase photoionization models. In 
this way, we can isolate the effects of the incident spectrum on species with ionization potentials in 
the range 1-5 Rydbergs (e.g., C III, C IV, Si III, Si IV). For example, Songaila & Cowie (1996a) 
argue that the observed increase in the Si IV/C IV ratio at z > 3.1 suggests that the He HI 
regions around ionizing sources may not have overlapped by that epoch, cutting off all radiation 
with hiz > 54.4 eV. Here, we primarily consider models for the averaged metagalactic background 
in which He H ionization fronts have overlapped before z ~ 4. We do consider a limiting case in 
which no photons above 4 Ryd are present in the background, as well as a case where we restore 
the much less heavily attenuated x-ray background. 

In the past decade, metal-line ratios in Lyman Limit cm“^) absorbers have been 

used to probe the shape of the ionizing spectrum, and to estimate the intrinsic metallicity of the 
absorbers (cf. Bergeron & Stasinska 1986; Sargent, Boksenberg &: Steidel 1988; Steidel, Sargent &, 
Boksenberg 1988; Steidel & Sargent 1989; Bergeron & Ikeuchi 1990; Miralda-Escude & Ostriker 
1990; Donahue & Shull 1991; Madau 1992). With the advent of new spectrographs on large 
telescopes, the lower Nhi absorbers can be exploited in the same way. These absorbers hold a 
special interest, because the source of their metal enrichment is still unclear. While the metal lines 
in Lyman limit absorbers have been associated with the haloes of bright galaxies, the number of 
Lya forest absorbers greatly exceeds the corresponding number of bright galaxies, and the Lya 
absorbers do not seem to show the same clustering properties as galaxies. 

As yet, metal lines have only been associated with “high-column” absorbers [N(H I) > 10^^'® 
cm“^], which numerically are a small fraction of the Lya forest. It is still possible that these 
clouds may represent an overlap of the high-column end of the pristine Lya forest clouds and the 
low-column end of the metal enriched clouds associated with galactic haloes. As a result some, 
if not all, of the metal-line absorbers may still be associated with the galaxies that contain the 
stars responsible for their metals. For example, over a Hubble time at z = 3, galactic outflows at 
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SOOVsoo km s ^ will transport heavy elements a distance 


d < (300 kpc)V3oo/i ^ 


1 + z 


- 3/2 


( 1 ) 


for a Friedmann universe with Hq = (100 km s~^)h and flo = 1- In practice, the period of heavy 
element injection may have lasted only 10® yrs, and the metals could move distances of only 30-50 
kpc from their sources (bright galaxies, dwarf galaxies, globular clusters). 

As Madau & Shull (1996) have emphasized, metal enrichment of the Lya absorbers implies a 
substantial Lyman continuum (LyC) emission accompanying the star formation at 2;^3.5. Thus, 
a large fraction of the metagalactic background at high redshift could be due to massive stars, 
a point also made by Giroux & Shapiro (1996) and Savaglio et al. (1997). Furthermore, the 
absorbers with metal lines may be sufficiently close to the local sources that the incident radiation 
field is dominated by local hot-star radiation. In contrast, if the metal enrichment of the Lya forest 
arises from Population III stars at z > 10 (Couchman 1985; Ostriker & Gnedin 1996), the stellar 
ionizing radiation will be greatly attenuated by intergalactic absorption, and the metagalactic 
background at z = 3 will be dominated by the harder spectrum of quasars. Because metal-line 
systems would then be less likely to be associated with nearby sources of stellar radiation, their 
absorption-line ratios would be more representative of the metagalactic ionizing background. 

Preliminary analyses (Songaila & Cowie 1996a; Savaglio et al. 1997; Haehnelt 1997) of the 
Lya forest clouds have assumed photoionization by simple power-law ionizing spectra. The only 
way these models account for filtering of the ionizing source spectrum by intervening clouds 
is to include a decrement at uhsII (4 Ryd). More realistic models include several additional 
effects. First, when absorption due to intervening clouds is included, the background spectrum 
just above the H I and He II ionization edges shows a flatter power law than that of the sources 
(Miralda-Escude & Ostriker 1990; Madau & Meiksin 1994; Giroux et al. 1995). Second, Haardt 
& Madau (1996) and Fardal, Giroux, &: Shull (1997) show, using cosmological radiative transfer, 
that cloud emission by recombination to the H I and He H ground states and He H Lya and 
2-photon radiation alters the extent of the decrement at I'Heii and changes the shape of the 
spectrum between 1 and 5 Rydbergs. Including a more realistic spectral shape will alter the 
relative populations of metal ions with edges in the 1 to 5 Ryd range even if the ionizing sources 
are all assumed to be quasars with intrinsic power-law ionizing spectra. 

In this paper, we consider the extent to which these alterations of the assumed shape of the 
incident radiation field affect the interpretation of the metallicity of these absorbers and the nature 
of the photoionizing sources. We use the data sets of Songaila & Cowie (1996a) and Savaglio et 
al. (1997), which include Si IV and C IV lines associated with higher column density Lya forest 
clouds. Songaila & Cowie (1996a) provide an important subsample for which C H lines are also 
observed. 
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2. The Metagalactic Radiation Background 
2.1. Ionizing Sources 

Assumptions about the nature of the ionizing sources of the metagalactic background are 
uncertain, since the intrinsic EUV spectra of quasars are not well dehned. For simplicity, we 
assume that all quasars have constant power-law spectra F^, oc Our choice of spectral index, 

as = 1.8, is consistent with HST observations of the EUV spectrum of quasars (Zheng et al. 1997) 
for radio-quiet quasars at 2 : = 1.5 — 2.0. As discussed by Fardal et al. (1997), if the ionizing 
background (1-5 Ryd) is dominated by these quasars, the He II Lya absorption optical depth 
THeii ~ 1 at 2 : = 3. For clarity, we denote the cloud-filtered, quasar-dominated mean intensity by 

JqjH- 

The characteristic spectrum of stellar sources in starburst galaxies is also not well known. 
Most models agree (e.g., Bruzual & Chariot 1993) that few photons are present with frequency 
above 54.4 eV, and that the spectrum falls off above hi/ ~ 45 eV. Recent line-blanketed, non-LTE 
hot-star models (Gabler et al. 1992; Sellmaier et al. 1996; Schaerer & de Koter 1997) vary 
considerably in the amount of radiation in the He I continuum (24.6 < hu < 54.4 eV) they predict, 
compared to the LTE line-blanketed models of Kurucz (1992). An ongoing uncertainty, for our 
purposes, is the lack of low-metallicity stellar models, which are more appropriate for high-redshift 
starburst galaxies. For the present work, we use a starburst spectrum (Sutherland & Shull 1997) 
discussed in the context of the ionizing spectrum of metal-forming galaxies by Madau & Shull 
(1996) [see their Fig. 3]. In brief, it is a time integrated spectrum of a Gaussian starburst of width 
2cj = 4 Myr containing 5000Mq of stars between 8 and 85 Mq with IMF slope T = 1.6 where 
dN{> M)/dM oc M~^ and N{> M) is the number of stars with mass greater than M. 

Estimates of the mean metagalactic ionizing intensity at 2 : > 2, based on the statistical 
diminution in the number of Lya forest lines close to quasars (the proximity effect), indicate a 
specific intensity at the hydrogen Lyman limit, uh^ of J_ 2 i = 0.7 — 1.5 where J_ 2 i is in units 
of 10“^^ergs cm“^ s“^ Hz“^ sr“^ (Gooke, Espey, & Garswell 1997). If all Lya absorbers possess 
a metallicity oi Z = O.OIZq by redshift 2 : ~ 3.5, the stellar ionizing photons associated with the 
production of the metals may be comparable to the ionizing photons produced by quasars (Madau 
& Shull 1996; Giroux & Shapiro 1996). If these ionizing photons are produced at higher redshift 
(e.g.. Population HI stars at 2 ; > 5) they will be highly attenuated by IGM absorption at 2 : > 4. 
A further uncertainty is the fraction, (/esc), of ionizing photons that successfully escape from the 
galaxy in which the stars reside (Dove & Shull 1994). It does seem possible, though, that stellar 
sources could be an important contributor to the metagalactic background, which could dominate 
in regions local to starburst galaxies. 

For example, the total Lyman continuum production rate of the Milky Way was recently 
inferred from GOBE observations of [N H] 205 /xm emission to be 3.5 x 10^^ photons s“^ (Bennett 
et al. 1994), and most starburst galaxies exceed the value S = ( 10^^)554 photons s“^. Assuming 
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that a substantial fraction, (/esc) > 0.1, of these photons escape these starburst galaxies, the mean 
photoionization rate of hydrogen due to the galactic source may exceed J_ 2 i = 1 if the absorber 
lies within a distance 


i?cr = (88 kpc) (/esc)^/' 


Us + 3] 
4.8 


( 2 ) 


Here, we have assumed the QSO spectra in the range 1-5 Ryd to be power laws with spectral 
index as ~ 1 . 8 . 


In this paper, we account for a stellar contribution to the mean intensity falling upon our 
absorbers in a global way. We assume that up to 2/3 of the ionizing emissivity from primary 
sources is stellar, as opposed to secondary sources due to recombination radiation from clouds. 
The remainder of the primary emissivity is then from AGN. We then compute the radiative 
transfer of primary and secondary sources through the IGM to derive the metagalactic ionizing 
background, JsQj{^)- The mean intensities Jqj (pure QSO sources) and JsQj (stars + QSOs) 
therefore bracket our estimates of the contribution of stellar sources to the uniform metagalactic 
background, Jmb{^)- However, even if the metagalactic ionizing background is quasar dominated, 
the absorber may primarily see a local unfiltered source of stellar photons, Jl(zz). We consider 
models which assume Jl/Jmb = 5, 10,20 at zz = bh for Jmb = Jqj and Jl/Jmb = 3,5, 10 , 20 for 
Jmb = JsQj- 


2.2. He II Ionization Fronts 


Within a picture of ionizing sources turning on in a neutral medium, Miralda-Escude & Rees 
(1994) made the important suggestion that the He H I-front which propagates away from sources 
with a sufficiently soft ionizing spectrum can lag the hydrogen I-front. Madau & Meiksin (1994) 
showed that, if as = 1.9, the overlap of He H I-fronts could be delayed to as late as z = 3. As a 
result, both Jakobsen et al. (1994) and Songaila & Gowie (1996) suggested that, at z = 3.1, He H 
I-fronts around the dominant ionizing sources may not have overlapped. Although a full treatment 
of overlapping He H I-fronts is beyond the scope of this paper, we can make a simple estimate of 
the ratio of I-front velocities for H I and He H. From the flux-limiting equations governing I-front 
propagation (Shapiro & Giroux 1987; Donahue &: Shull 1987), the velocity relative to the local 
Hubble flow is given by 



Si 

ATrrfni ’ 


( 3 ) 


where i refers either to the H I ionizing front (1 Ryd continuum) or to the He H ionizing front 
(4 Ryd continuum). (The He I ionizing front is assumed to be coincident with the H I front, 
and we neglect recombinations and attenuation within the ionized zone.) Here, Si (photons s“^) 
represents the photon production rates of the source in the H I and He H continua, n* is the 
density of H I or He H, and Vi is the front distance measured from the source. 



- 6 - 


The propagation speeds of the H I and He II ionization fronts will be governed by the 
“flux-to-density” ratio, (Si/riirf), at 1 and 4 Ryd. For a QSO spectrum with power-law index a*, 
one finds that <S'Heii/'S'Hi = 4“"". The primordial helium abundance by mass has been estimated at 
values Yp = 0.231 ± 0.006 (Skillman & Kennicutt 1993) and Yp = 0.232 ± 0.003 (Olive & Steigman 
1995), while Copi et al. (1994) have suggested a 2cj concordance range of Yp = 0.221 — 0.243 
from theoretical models of Big-Bang nucleosynthesis. If helium has a cosmological abundance, 
nue/nu = 0.0785 by number (Y = 0.239 by mass), the flux-to-density ratios at a fixed radius r* 
are equal for 


a, = —- 


ln[ raHe/nH] 
In 4 


1.84 


( 4 ) 


At this critical spectral index, the H I and He H ionization fronts will coincide. Interestingly, this 
critical index is very close to the mean index found by Zheng et al. (1997) for radio-quiet quasars. 
For AGN with harder spectra (a^ < 1.84), the He H ionization front will precede the H I front by 
a small amount, equivalent to a few optical depths in the 4 Ryd continuum. However, for those 
AGN with softer spectra (a^ > 1-84), the He H front will propagate at a lower speed and lie well 
within the H I front. For example, Zheng et al. (1997) found that radio-loud quasars had a mean 
index ag ~ 2.2. In this case, the regions outside the He H fronts will see no 4 Ryd continuum 
radiation. 


We therefore consider a limiting case in which all photons above 4 Ryd are removed from the 
metagalactic background spectrum. We also consider an intermediate case, arising from the fact 
that quasars emit radiation at least into the x-ray range of energies, so that high-energy photons 
will not be strongly attenuated, even in regions where helium is entirely in the He H ionization 
stage. In the Lya clouds, we estimate the effect of the presence of these high-energy photons by 
presuming that H H regions have overlapped by z = 4, and that the Lya forest clouds, whose 
distribution in A^hi is known (see §2.2), have a maximal He H/H I ratio. In interpreting He H 
Gunn-Peterson measurements, several groups (Madau & Meiksin 1994; Giroux et al. 1995) find 
that the ratio, ry = N(He H)/N(H I), must be 50-100 in order to explain the A < 304 A absorption 
as He H line opacity in the H I Lya forest. This large ratio arises because He H is more difficult 
to photoionize than H I. 

The H neutral fraction fpi is always less than (nHe/f^n)??”^) where ry = N(He H)/N(H I) 
within an absorbing cloud. Thus, the maximal ratio rjmax is bounded by limits on the baryon 
density from Big Bang nucleosynthesis (BBN), < 0.024 (Gopi et al. 1994). At z = 3.4, if 
only absorbers with Npi < 10^® cm“^ are considered to contribute to we find r]raax < 5000. 
This limit produces a local continuum optical depth at 4 Ryd, dr/dz > 180 at z = 3.4, falling to 
Ahy V = 25z//f, the same continuum optical depth as that at 1 Ryd. As a result, for this model 
with non-overlapping He H I-fronts, we only remove photons between 4 and 25 Ryd, retaining 
metagalactic radiation background for v > 25 i'h- This is a more conservative estimate than if all 
He is assumed to be in He H and distributed uniformly. 
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2.3. Cloud Opacity/Emission 


For a random spatial distribution, the local continuum optical depth of the intervening 
absorbers is given by (Paresce, McKee, Sz Bowyer 1980) 


dr(zv) _ d‘^Nc 

dz Jo dNnidz 

where, for clouds composed of H and He, 

Nuel 


[1 - dNni , 




INhi\ 


CTHeli^) 


NhsII 

[ Nhi J 




( 5 ) 


( 6 ) 


and is the distribution of absorbers in column density and redshift. This distribution has 

often been parameterized by a function of the form 


dNffjdz 


= A{l + zrN^^ 


( 7 ) 


for Ni{}l I) < NH[{cm~‘^) < Nu{}i I). Using new line lists from Keck data, Fardal et al. (1997) 
have completed a review of the statistics available on Lya absorbers, Lyman Limit Systems, and 
damped Lya systems. They suggest that the distribution of absorbers is best represented by a 
broken power law in Nhi- If 


dNndz 


Ai{l-\- z ^, 


( 8 ) 


where A^i 7 = Nhi/{^^^'^ cm“^), 7 = 2.585, and Pi and Ai represent coefficients appropriate for a 
given range in A'"i 7 (see Table 2). We adopt their Model 2 for this paper. Use of either of their 
other two models does not affect our conclusions. 


Just as absorption due to intervening clouds will strongly alter the background radiation 
spectrum, the re-emission of ionizing radiation from the clouds will also affect the radiation 
(Haardt & Madau 1996; Fardal et al. 1997). Recombinations in the clouds produce continuum 
radiation (H I Lyc, He H Lyc, He H Bac, and He H 2-photon) as well as He H Lya (A304) line 
radiation. Accounting properly for the Lyc and Bac radiation is critical to the determination of the 
level of the mean intensity above 4 Ryd. While the line and 2-photon radiation are not important 
for that, they strongly alter the spectrum between 1 and 3 Ryd and affect the populations of ions 
with thresholds in that range (see Table 1). A complete discussion of the treatment of reemitted 
radiation used in computing the filtered spectra is given in Fardal et al. (1997). 


3. Lya Forest Cloud Models 

As stated previously, we neglect the possibility that the absorbers consist of more than one 
thermal and ionization phase to concentrate on the effects of the shape of the photoionizing 
radiation spectrum. We model the Lya forest absorbers using the photoionization code CLOUDY 



version 90 (Ferland 1996). The absorber clouds are modelled as plane-parallel slabs of constant 
density, illuminated on both sides by the ionizing radiation field. We adopt a fiducial column 
density N(H I) = 10^® cm“^ and assume a metallicity Z = 10’ -2^0. We adopt the spectral 
shapes for the incident radiation discussed in §2.1 and vary the strength of the incident radiation 
through the ionization parameter, U = n-fjnH, where and nn are the number densities of 
ionizing photons and hydrogen respectively. Each curve in Figs. 2, 4, and 6 represents a set 
of photoionization models paired with the corresponding spectrum shown in Figs. 1, 3, and 5, 
respectively. 

Although we focus this paper on the effect of the spectral shape, we have also explored the 
effect of varying the parameters assumed above. For the spectral shapes considered here, there 
is little change in the ratios (C II/C IV is reduced by ~ 10% and Si IV/C IV is raised ~ 10%) 
if N(H I) is increased by a factor of 10, as long as the temperatures for the clouds are assumed 
to be consistent with the thermal equilibrium solution of the CLOUDY model. An increase (or 
decrease) by a factor of 3 in the assumed metallicity has a slightly larger effect. The C II/C IV 
fraction is reduced 20 — 30% with a factor of 3 reduction in metallicity and raised 8 — 15% by a 
factor of 3 increase in metallicity. The ratio Si IV/C IV is raised as much as 10% and reduced 
by 6% or less with the same changes in metallicity. Constraints on the H I columns of the clouds 
preclude large reductions in the assumed density for one-phase photoionization models, to avoid 
excessive length scales for the clouds and to ensure that contributions to the baryon density from 
these N(H I) ^10^^ cm“2 absorbers do not exceed BBN limits. 

We normally allow the temperature to be that solved for by CLOUDY, assuming thermal 
equilibrium between photoelectric heating and radiative cooling. If the clouds are actually cooler, 
for example if they have expended some energy in expansion, or if they retain thermal memory 
of earlier epochs when Compton cooling off the microwave background was important, the ratios 
Si IV/C IV vs. C II/C IV may be increased by as much as a factor of 2 at low values of C II/C IV 
(high values of U). If, as suggested by Haehnelt, Rauch, & Steinmetz (1997), the temperatures 
of the clouds are larger than that expected from photoionization thermal equilibrium, the Si IV 
fraction is more temperature-sensitive than the C IV fraction, and Si IV/C IV is decreased further 
relative to C II/C IV. These higher temperatures, suggested as a way to increase the apparent 
thickness of the clouds, also imply much smaller metallicities, as the C IV fraction is decreased 
much less compared to the H I fraction by higher temperatures. For example, adopting the Jqj 
spectrum and assuming U = 10“2^ -we find that increasing the cloud temperature from 25,000 K to 
70,000 K decreases the C IV fraction by less than 20% but reduces the H I fraction (and inferred 
metallicity) by a factor of 8. 

There is a further effect if the assumed temperature of the cloud is increased to the point that 
collisional ionization of H I becomes important. In this case. He II/H I may be larger than that 
assuming pure photoionization only, N(He II) may exceed l/aneii-, and the cloud may become 
self-shielding. Even if this is not the case for a cloud with N(H I) = 10^^ cm“2^ this becomes 
increasingly likely for greater values of N(H I), so that radiative transfer within the individual 



- 9 - 


clouds becomes important. In Fig. 7 we show the results of models that explore some of these 
temperature effects but defer a complete discussion to a later paper. 


4. Results 

4.1. The Metagalactic Radiation Background 

Figure 2 summarizes the results of our photoionization models which assume different sources 
of photoionizing radiation for the metagalactic background. The solid and dotted curves both 
assume photoionization by sources that possess the same power-law spectral shape (a^ = 1.8) 
but differ in the treatment of the filtering of their radiation by intervening clouds (see Fig. 1). 
Models incorporating a pure-AGN spectrum, Jqj, indicate that the Si IV/C IV ratio may be 
enhanced by almost a factor of 2 over the simpler model for the absorbers with C II/C IV ^0.1. 
Even in this regime, Si is still overabundant relative to C, as previously argued by Songaila & 
Cowie (1996a). However, a factor of two enhancement may be sufficient. If the absorbers are 
photoionized, an increasing C II/C IV ratio is directly related to a decreasing ionization parameter 
(Figure 2). At high C II/C IV (low U) there are few high-energy photons available to ionize the 
higher ionization stages of silicon. As a result, the enhanced ionization of Si III by He H Lya 
photons emitted by clouds is enough to increase the ratio of Si IV/C IV over that of the broken 
power-law models. Compared to more realistic calculations with cloud re-emission, models for the 
spectral shape of the metagalactic background with a simple break at 4 Ryd underestimate the 
number of high-energy photons in the metagalactic background. At low C II/C IV (high U), the 
Si IV fraction is strongly depleted by the increased amount of high-energy photons in the realistic 
spectrum. 

A metagalactic ionizing background dominated by stellar sources, Jmb = JsQj-, increases the 
Si IV/C IV ratios at C II/C IV > 0.1 by 20 — 50% over that of the Jmb = Jqj case. Once again, 
using a more realistic spectrum increases the Si IV/C IV ratios by almost a factor of 2 compared 
to spectra with simple attenuation of He H ionizing photons. As in the pure ACN case, the simple 
representation overestimates the corresponding intensity at higher energies. To be compatible with 
observed Si IV/C IV ratios if C II/C IV > 0.1, the relative abundance of Si must be increased by 
approximately a factor of 2. In all cases, the large Si IV/C IV ratios in absorbers with C II/C IV 
<0.1 are difficult to explain in the context of single-phase models unless Si/C is enhanced by an 
order of magnitude. Chemical evolution models with Si/C > 3 (Si/C)© are unrealistic, even if the 
nucleosynthesis is dominated by massive stars (Woosley &: Weaver 1995). 
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4.2. A Photoionizing Background Dominated by Local Sources 

One resolution to the puzzle of large Si IV/C IV at small C II/C IV may be that the ionizing 
radiation incident on the absorbers is dominated by stellar sources of radiation. As Fig. 4 shows, 
if the ionizing spectrum is entirely stellar, very large ratios of Si IV/C IV are always possible 
due to the steep dropoff in radiation more energetic than 45 eV. In practice, however, as Figs. 

3 and 4 show, if the metagalactic background is primarily due to quasars {Jmb = Jqj)i even if 
Jl/Jmb = 20, high Si IV/C IV ratios at low C II/C IV require Si enhancements of a factor of 10. 

If the metagalactic background is dominated by stellar sources {Jmb = JsQj)-, large 
enhancements of Si may be unnecessary (see Figs. 5 and 6). For many absorbers with C II/C IV 
> 0.1, no overabundance of Si is necessary if Jl/Jmb = 3 — 10. At low C II/C IV, Jl/Jmb 
must exceed 10 if enhancements in Si exceeding 2 are to be avoided. As our order of magnitude 
estimates in §2.1 indicate, this would require that the absorbers lie within 30 kpc of starburst 
galaxies. However, this situation is quite feasible, given realistic constraints on absorber size and 
metal transport distances. 


4.3. Non-Equilibrium Temperatures and Non-Uniform Radiation Fields 

As Figs. 7a and 7b show, if the temperatures of the clouds are cooler than the photoionization 
thermal equilibrium temperatures (Ferrara & Giallongo 1997; Zhang et al. 1997), much higher 
ratios of Si IV/C IV are possible for low C II/C IV (high U) clouds. These high U clouds 
are precisely the clouds for which the photoionization thermal equilibrium temperature is 
highest (T w 40,000 K). The models denoted by the dotted curves in Figs. 7a and 7b assume 
T = 15,000 K, which may not be compatible with the observed line widths, and which are cooler 
than is usually assumed for Lya forest clouds. If higher temperatures are associated with the 
clouds, for example the dashed curves in Figs. 7a and 7b, which assume T = 50,000 K, it is more 
difficult to account for the ratios in a one-phase model. 

Another way to enhance Si IV/C IV at low C II/C IV, as suggested by Songaila & Cowie 
(1996a) and Savaglio et al. (1997), is to assume that no photons above 4 Ryd are present in the 
incident ionizing spectrum. This is the limiting case of a cloud embedded in a region where He 
H I-fronts have not yet overlapped. The Si IV/C IV ratios in Fig. 8 are not as high as those for 
the limiting case when all incident radiation is starburst radiation, since radiation between 45 eV 
and 54 eV is included, which will ionize Si IV to Si V. This effect is seen in comparing models 
which assume pure AGN and starburst/AGN sources in Fig. 8. As Fig. 8 also shows, a careful 
treatment of the propagation of these He H I-fronts is necessary, particularly because most AGN 
possess x-ray emission with spectra that flatten to approximately above 0.3 keV. We have 
shown that including high-energy radiation with E > 25 Ryd is sufficient to reduce Si IV/C IV 
ratios substantially at low C II/C IV. 
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5. Discussion 

From the results of our one-phase photoionization models we draw the following conclusions: 

1. Ionizing Radiation Field- The metagalactic radiation field is likely to include both 
power-law (AGN) and stellar (hot-star) components. The Lya forest absorbers with metal lines 
may also experience a local radiation field from starburst galaxies within 50-100 kpc. 

2. Si/C Overabundance - For plausible mixtures of stellar and AGN spectra in the 
metagalactic background, our photoionization models produce enhanced Si IV/C IV ratios, 
consistent with high-z absorbers with Si/C ~ 2(Si/C)Q at low ionization (C II/C IV > 0.1). For 
absorbers with C II/C IV < 0.05, it is difficult to account for the high values of Si IV/C IV unless 
Si/C > lO(Si/C) 0 , an unrealistically large value for massive-star nucleosynthesis. These systems 
may include photoionization from local stellar sources as well as hot, collisionally ionized gas. 

3. Local Ionizing Sources - If the radiation field incident on the absorbers is dominated by a 
nearby starburst galaxy, the Si IV/C IV ratios are further enhanced. If C II/C IV > 0.1, no Si/C 
overabundance is necessary to explain Si IV/C IV ratios if absorbers are within about 40 kpc of 
a starburst galaxy and the background is dominated by stellar sources. Even if the metagalactic 
background is dominated by AGNs, close proximity to a starburst galaxy may reduce the needed 
Si/C overabundance to a factor 1.5. When C II/C IV ^0.05, it remains very difficult to account for 
the high values of Si IV/C IV with photoionized models, although models with locally dominated 
radiation may only require Si/C enhancements of a factor of 2-3. 

4. Non-Overlapping He II I-Fronts - Around quasars whose ionizing continua have steep 
spectral indices (a* > 1.84), the He II I-fronts will lie within the H I I-fronts. If many absorbers 
lie in regions where all photons above 4 Ryd are attenuated (Songaila & Cowie 1996a), we obtain 
good agreement with almost all measurements of Si IV/C IV if Si/C is enhanced by a factor of 
2-3. However, including even a small contribution of higher energy photons in the background 
increases the needed Si/C overabundance to an order of magnitude for absorbers with the lowest 
C II/C IV ratios. 

5. Temperature Effects - If the absorbers are cooler than expected for thermal equilibrium 
between photoelectric heating and radiative cooling, the Si IV/C IV ratio is increased for a 
given C II/C IV ratio. At low C II/C IV, this enhances Si IV/C IV by a factor of 5 in our 
photoionization models. Conversely, if the absorbers are hotter, Si IV/C IV is lower for a given 
C II/C IV ratio. 


While we have described several processes that might increase the Si IV/C IV ratio, 
an important additional constraint on these possibilities is the fact that the high ratios 
are preferentially found at z > 3.1. Recent observations (Boksenberg 1997) challenge this 
interpretation, by finding high Si IV/C IV ratios in absorbers at z = 2 — 3. If, however. Si IV/C IV 
rises above z > 3.1, this argues for a time dependence to whatever process increases this ratio. It 
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is this property, as well as the increased He II absorption toward Q0302-003 at z = 3.28 (Jakobsen 
et al. 1994), that Songaila & Cowie (1996a) use to support their suggestion that the absorbers 
lie in regions where He H I-fronts have not overlapped at z > 3.1. From a higher resolution 
HST/GHRS spectrum of Q0302-003, Hogan, Anderson, & Rugers (1997) find evidence for residual 
transmitted flux below the He H edge. Their 95% confidence upper limit, thbII < 3, makes it less 
necessary to propose that Q0302-003 lies in a region where He H I-fronts have not overlapped, but 
does not preclude the suggestion that such regions existed at z > 3.1. If this is the case, radiation 
outside of He HI regions would largely limit the range in ionization stages for carbon to H-IV, and 
for silicon to H-V. In AGN, however, the large fluxes of unattenuated photons above 45 eV can 
ionize Si IV and C HI, so that high observed Si IV/C IV ratios would be difficult to explain for 
large U (small C H/C IV). As a result, the spectrum below 4 Ryd must be dominated by stellar 
sources. This may be difficult to achieve with the known quasar luminosity functions at z > 3.5. 
In addition, soft x-ray radiation from quasars, which is less attenuated by intervening clouds even 
in the low He HI porosity case, may again make high Si IV/C IV, low C H/C IV absorbers difficult 
to understand if they are solely photoionized by a metagalactic background. 

New or improved measurements of He H absorption may soon indicate whether He H I-fronts 
have overlapped well before z = 3.1. In that case, a sharp change in the shape of the metagalactic 
radiation field may be less plausible. Still, in general, the trend is likely to be more attenuation of 
photons above the 4 Ryd limit with increasing redshift. Other time dependent effects that enhance 
the ratio of Si IV/C IV are possible. Higher abundances of Si relative to C are associated with 
metal yields from the most massive stars. However, in a flat universe with h = 0.75, the age of the 
universe exceeds 10® years by z = 3.1, time enough for lower mass stars to enrich the gas with 
carbon. For example, if multiple supernovae eject metal-enriched gas into the IGM, the additional 
4 X 10® years between z = 3.5 and z = 2.5 may make it more likely that the gas surrounding the 
high mass stars has been enriched with carbon from a previous episode of star formation. 

The temperature effects we discuss in §4 may also have a time dependence. If the IGM has 
been uniformly enriched by a much earlier episode of Population HI star formation, and if the 
Lya forest clouds with metals have their origin in growing overdensities in the IGM, then the 
absorbers observed at higher 2 : may have formed at an earlier epoch. If Compton cooling off the 
cosmic microwave background was a dominant coolant at this epoch (Miralda-Escude & Rees 
1994), the cloud may have retained memory of this lower temperature. This effect is probably not 
large, since there is not a significant evolution in the observed line widths. 

Alternatively, one may resolve the puzzle of the high Si IV/C IV, low C H/C IV absorbers 
by relaxing the single-phase model which we used to explore the effects of radiative transfer. As 
we have emphasized in a previous paper (Giroux et al. 1994), the production of heavy elements 
in QSO absorption systems is naturally accompanied by hot gas due to supernovae and hot-star 
winds. 


We thank Mark Fardal for the calculated metagalactic background radiation spectra. We 
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Table 1. Significant Multi-stage Metals 


Species ^th(Ryd) 


C II, 

III, IV 

24.38, 

47.89, 

64.49 


NI, 

II, III, IV 

14.53, 

29.60, 

47.45, 

77.47 

0 1 , 

II, III 

13.62, 

35.12, 

54.94 


Ne I, 

II, III 

21.57, 

40.96, 

63.46 


Si IL 

, III, IV, V 

16.35, 

33.49, 

45.14, 

166.77 

A1 II 

,III 

18.83, 

28.45 



S II, 

III, IV, V 

23.33, 

34.83, 

47.3, 

72.68 

Fe II 

, III, IV 

16.19, 

30.65, 

54.8 



Table 2. Absorber Distribution 


Range in V 17 

A, 


10“® - 10“3 

0.132 

1.41 

10“3 - 0.02 

6.31 X 10“® 

1.85 

0.02 - 10® 

2.9 X 10“2 

1.46 
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Fig. 1.— Mean intensity, in relative units, versus energy in Ryd. All curves have been normalized 
to the same values at = uh (1 Ryd). Solid curve assumes IGM-filtered sources with AGN 
(j^- 1 . 8 ) spectrum at z = 3.4. Dotted curve assumes a mean intensity with above 

(1 Ryd) and v = VHeii (4 Ryd), with decrease by a factor of 2 above i^Heii, designed to give 
same ratio, J(4 Ryd)/J(lRyd), as the properly filtered case. Short-dashed curve assumes an IGM- 
filtered spectrum; 2/3 starburst (from Sutherland & Shull 1997) and 1/3 AGN . Long-dashed curve 
assumes above u = vh and u = izHeii, and a decrease by a factor of 12 above i^Heii, to give 

same ratio, J(4 Ryd)/J(lRyd), as the mixed-source spectrum above. 

Fig. 2.— Ratios Si IV/G IV versus C II/C IV in photoionized models with N(H I)= 10^^ cm“^. 
Curves match the spectra of Fig. 1. Triangles are from measured column densities of Si IV, C IV, 
and C II (Songaila & Cowie 1996a,b). Squares are measurements and upper limits from Songaila & 
Cowie (1996a,b). Diamonds are measurements and upper limits from Savaglio et al. (1996). Same 
data are also used in Figs. 4, 6, 7, 8 with model calculations of Si IV/C IV versus C II/C IV. 

Fig. 3.— Mean intensity, in relative units, versus energy in Ryd. All curves have been normalized 
to the same values at = vh- All curves are mixtures where Jmb is the IGM-filtered AGN 
spectrum shown Figure 1, and Jl is the unfiltered (“local”) stellar spectrum (Sutherland & Shull 
1997). Curves assume Jl/Jmb = 0 (solid), 5 (dotted), 10 (short-dashed), and 20 (long-dashed). 
Dot-long-dashed curve assumes only local stellar radiation, Jl. 

Fig. 4.— Ratios of Si IV/C IV versus C II/C IV in photoionized models with N(H I)= 10^® cm“^. 
Curves match the spectra of Fig. 3. 

Fig. 5.— Mean intensity, in relative units, versus energy in Ryd, where all curves have been 
normalized to the same values at v = vh- All curves are mixtures where Jmb is the filtered 
stellar/AGN spectrum shown in Figure 1, and Jl is the unfiltered stellar spectrum from Sutherland 
& Shull (1997). Curves assume Jl/Jmb = 0 (solid), 3 (dotted), 5 (short-dashed), 10 (long-dashed), 
and 20 (dot-dashed). 

Fig. 6.— Ratios Si IV/C IV versus C II/C IV in photoionized models with N(H I)= 10^^ cm“^. 
Curves match the spectra of Fig. 5. Dot-long-dashed curve assumes only local stellar radiation, 
Jl- 

Fig. 7.— (a) Effects of varying cloud temperature. All cloud models assume a metagalactic 

background with quasar sources only, Jmb = JQ,f- Solid curve adopts the equilibrium temperature 
computed by CLOUDY for a given ionization parameter and spectral shape (varies between 
T = 16,500 K and T = 45,000 K.) Dotted curve assumes constant T = 15,000 K for all models; 
dashed curve assumes constant T = 50,000 K for all models, (b) Same as 7a, except background 
includes both quasars and hot stars, Jmb = JsQj- 

Fig. 8.— Ratios Si IV/C IV versus C II/C IV in an IGM in which He II ionization fronts have not 
overlapped. Solid curve repeats models for the filtered starburst/AGN spectrum in Figure 2. Dotted 
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assumes filtered starburst/AGN spectrum with all photons with E > A Ryd removed. Long-dashed 
curve assumes AGN spectrum with all photons with > 4 Ryd removed. Short-dashed curve 
assumes hltered starburst/AGN spectrum with 4 < i? < 25 Ryd photons removed. 





Si IV/C IV 




Log J 



1 


10 

Energy (Ryd) 



Si IV/C IV 






Si IV/C IV 




Si IV/C IV Si IV/C IV 




Si IV/C IV 




